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ABSTRACT 

Practical techniques for designing transfer trajectories 
between Libration Point Orbits (LPOs) are presented. 
Motivation for development of these techniques was 
provided by a hardware contingency experienced by the 
Solar Heliospheric Observatory (SOHO), a joint 
mission of the European Space Agency (ESA) and the 
National Aeronautics and Space Administration 
(NASA) orbiting the LI point of the Sun-Earth system. 
A potential solution to the problem involved a transfer 
from SOHO’s periodic halo orbit to a new LPO of 
substantially different dimensions. Assuming the 
SOHO halo orbit as the departure orbit, several practical 
LPO transfer techniques were developed to obtain new 
Lissajous or periodic halo orbits that satisfy mission 
requirements and constraints. While not implemented 
for the SOHO mission, practical LPO transfer 
techniques were devised that are generally applicable to 
current and future LPO missions. 

1. INTRODUCTION 

To date there have been several spacecraft orbiting 
either the LI and/or L2 Lagrange, or libration, points of 
the Sun-Earth/Moon barycenter system. In the future 
we may expect many more solar and astrophysics 
missions to take advantage of the observational benefits 
associated with these locations. Renewed interest in the 
possible use of the coll inear points of the Earth-Moon 
system has arisen recently as well. The use of Libration 
Point Orbits (LPO) has been accompanied by a rich and 
growing literature on the design, modeling, and analysis 
of a variety of orbits in the restricted three-body 
problem (RTBP). Less attention has been paid, 
however, to the problem of transfers between distinctly 
different LPOs; in particular, the problem of transferring 
between two LPOs that do not intersect. The problem 
of how to design such transfers recently arose as a 
practical problem for a currently operational LI 
mission. The question was posed: how does one design 
such a transfer, and conduct it efficiently with the 
existing operational support software and with regard to 
mission constraints including remaining fuel? 

The event spurring this question was a contingency 
occurring with the ESA/NASA Solar and Heliospheric 
Observatory (SOHO), history’s second and longest 
serving LI orbiter (launched December 1995). In May 
2003, trouble with the High Gain Antenna (HGA) 


gimbal motor led to the HGA being left in a fixed 
position. This condition made high data rate contacts 
unachievable over long portions of a revolution, which 
jeopardized not only the science return but the 
continued viability of the mission itself. In response, a 
study was conducted to investigate the possibilities for 
transferring SOHO to an LPO of different dimensions in 
hopes of improving the HGA-to-Earth coverage 
geometry. 

The investigation took as its starting point the 
established SOHO LPO — the periodic type commonly 
referred to as a ‘halo’ orbit — and considered transfers to 
both halo-type and the non-periodic Lissajous-type 
orbits of distinctly different dimensions. Though 
techniques are possible involving a return to the Earth 
for a lunar swingby-assisted transfer to a different LPO, 
the work described here is limited to transfers conducted 
entirely within the LI region. Practical transfer solutions 
were constructed for both types and found to be 
affordable fuel-wise but not advisable for SOHO due to 
past anomalies that have seriously degraded its 
capability to perform orbit maneuvers of large size [1]. 
Nevertheless, practical LPO transfer techniques were 
devised that have general applicability. 

2. THE SOHO MISSION HALO ORBIT 

SOHO’s halo orbit is referred to as a Class 2 orbit, 
indicating that its sense of revolution about LI is 
counter-clockwise as seen from Earth. Both halo and 
Lissajous orbits are described with respect to the usual 
synodic reference frame of the RTBP, commonly called 
the Rotating Libration Point (RLP) frame in work at 
Goddard Space Flight Center (GSFC). The x-axis of the 
LI -centered RLP frame points toward the Earth-Moon 
barycenter, the r-axis points up toward the North 
Ecliptic Pole (NEP), and the ^y-axis completes the right- 
handed triad. The RLP system should be assumed 
herein, except for impulse variables for which the 
reference frame is described in Section 3. 

Because of solar radio interference constraints, SOHO is 
required to avoid the region within 4.5 degrees of the 
Sun called the Solar Exclusion Zone (SEZ). The 
minimum halo r-axis amplitude, A z , needed to satisfy 
the above SEZ requirement was 120,000 km, the value 
selected during mission design. Richardson shows in 



his development of the formulation for periodic halo 
orbit motion about the collinear points [2] that the non- 
linearities of the problem force a relationship between 
the r-amplitude and ^-amplitude. Richardson provides 
an algebraic expression that yields a ^-amplitude (Ay) of 
666,672 km when A z is 120,000 km. Also, because the 
motion in the xy plane is coupled, the corresponding jc- 
amplitude (A x ) for this case is 206,448 km. 

The nominal SOHO halo orbit is depicted in the RLP 
xy- and ^-projections in Fig. 1 and Fig. 2, respectively. 
For contrast, Fig. 1 and 2 also show the small-amplitude 
Lissajous orbit belonging to NASA’s Advanced 
Composition Explorer (ACE) mission. Additional 
information on the SOHO mission can be found in [1]. 



Fig. 1. SOHO Halo RLP XY Plane Projection 



Fig. 2. SOHO Halo RLP XZ Projection 


3. LPO TRANSFER TECHNIQUE 

Peculiarities of the SOHO anomaly led to two lines of 
inquiry, both of which were well suited to an initial 
exploration of the transfer problem. The first asked what 
is required to change Ay (in particular, to reduce Ay) 
while leaving A z unchanged? The second asked, what is 


required to change A z (in particular, to increase A z ) with 
no change required to Ay? 

The transfer problem is seen to have two main parts. 
First, a target LPO must be selected, which means that 
mission requirements, goals, and constraint criteria must 
be considered to arrive at an appropriate choice of target 
LPO type and overall dimensions. The second part of 
the transfer problem then consists of designing a 
transfer trajectory arc connecting the initial and final 
LPOs. The design of the final LPO also requires 
selection of a target point, or patch point, for the transfer 
arc to intersect. 

For the problem of target — or destination — LPO 
selection, two things are apparent at once. First, the 
target LPO should be of the same class as the initial 
LPO. Attempting to jump onto an LPO of the opposite 
class — that is, one with motion in the opposite 
direction — would necessarily be prohibitively costly. 
Second, a minimum of two impulses — an impulse to 
initiate departure from the original LPO and an insertion 
impulse on the target LPO — would seem necessary. 
Thus, in addition to selecting target LPO dimensions, 
departure and insertion locations must be chosen. 
Finally, coordinates (position and velocity components) 
for the target insertion location must be calculated. 

The trajectory propagation and targeting software used 
for this work is the same as that used for flight 
operations — namely, the program called 3 4 Swingby ’ [3]. 
Targeting within Swingby employs an iterative 
differential corrections (DC) process. The user 
constructs a targeting scheme consisting of independent 
targeting variables (components of the impulses) and 
dependent variables (targeting goals). Generally a 
scheme is constructed such that the number of 
independent variables equals the number of dependent 
variables to provide the differential corrector with a 
‘square’ problem to solve. Dependent variables may be 
Cartesian position components, velocity components, or 
a mixture. Independent and dependent variables both 
may be separated in space and time. Swingby 
propagates trajectories numerically with full operations- 
level force modeling invoked. 

Four trial cases were developed for both the Ay 
reduction and the A z amplification transfer problems. In 
all cases the impulse variables used are with respect to a 
spacecraft-centered delta-V coordinate frame. In this 
frame, the *-axis points to the Sun, the r-axis points 
toward the NEP, and the j-axis completes the right- 
handed system. All references to impulses (i.e., delta- 
Vs, or burns) herein are with respect to this system, 
which was originally defined for the SOHO mission. (It 
was naturally suited to SOHO’s solar-pointing 
orientation and thruster configuration, but is well suited 
to LPO work generally because its axes are virtually 
parallel with the RLP and related synodic systems.) 




3.1 A Y Reduction Transfer Design 

The question of reducing the Ay of a halo orbit is 
examined first. In the case of SOHO, its halo is very 
near the minimum Ay halo. Periodic halo orbits must 
have an A y > 654,276 km [4]. Smaller ^-amplitude 
LPOs are necessarily of the Lissajous type. Therefore, 
for SOHO, any substantial reduction of Ay necessarily 
means the transfer would be to a Lissajous orbit. 

In his development of the third-order solution to the 
equations of motion for LPOs, Richardson [2, 5] shows 
that for Lissajous orbits both A z and Ay are free 
parameters (unlike for periodic halo orbits). (For both 
halo and Lissajous orbits, A x and Ay are related via a 
proportionality constant.) In fact, to fully specify a 
Lissajous orbit, four parameters are required. In 
addition to A z and Ay, two phase angles (j> xy and 
y/f — pertaining to the in-plane and out-of-plane 
oscillations — must be specified. In short, specifying the 
amplitudes establishes the overall dimensions, while the 
phase angles establish a given point on the orbit. The 
selection of (j> xy and y/ z also establishes the ‘class’ of 
orbit, i.e., its sense of revolution about LI as seen from 
Earth. (The angle fay is measured in the RLP xy-plane 
clockwise from the -x-axis, and the angle y/ z is the 
phase with respect to the z-axis.) In summary, given the 
four parameters of the desired Lissajous, the third-order 
theory yields the state coordinates of the desired target 
insertion point. 

How to choose the target insertion point and the transfer 
arc duration become the next two questions. An 
immediate answer to the first was lent by the past 
practices of targeting both Earth-to-Ll transfers and 
LPO stationkeeping maneuvers to target coordinates on 
thejy = 0, or xr, plane. Swing by actually facilitates such 
targeting via a user-selectable propagator stopping 
condition feature that identifies RLP plane crossings. 
Experience with LPO stationkeeping and Lissajous 
control helped suggest an answer to the second. It 
seemed a minimum of one-half a revolution about LI 
would be necessary to keep transfer costs reasonable. 
The total period of revolution in the halo orbit is about 
178 days, though interestingly the flight time from the 
Earth-side xr-plane to the Sun-side xr-plane is 90.5 days 
while the return to Earth-side takes only 87.5 days. 

At the time this work began, SOHO was coming up to a 
crossing of the Earth-side xr-plane occurring on July 7, 
2003. It was this xr-plane impact point state that was 
used as the nominal halo departure point for all transfer 
cases developed. Thus, the approach devised is that a 
departure impulse is targeted on one or more coordinate 
goals (e.g., the position on the target LPO) at a specified 
future xr-plane crossing. It was decided to target three 
of the trial cases on the Sun-side xr-plane, i.e., Vz 
revolution from the Earth-side departure point. The 
fourth case targeted a return to the Earth-side crossing, 


making for a full revolution transfer arc. Once the xr- 
plane impact point is successfully achieved, the 
remaining step is to apply the LPO insertion impulse. 
The insertion impulse supplies the difference between 
transfer arc arrival velocity vector components and the 
three components of the target LPO velocity. 

Detailed development of the Ay cases is now considered. 
The first three cases — AY1, AY2, and AY3 — involve 
reducing Ay to 500,000 km (resulting in a Lissajous) 
using three targeting scheme variations. Since it was 
not the intent to change A z , the amplitudes used to 
construct the Lissajous were A z = 120,000 km and Ay = 
500,000 km. Cases AY1 and AY2 involved Vz 
revolution transfers to the Sun-side xr-plane crossing, 
hence to complete specification of the target Lissajous 
insertion state the selected phase angles were fay = 0 and 
y/ z = -90 (herein, units for the phase angles are 
understood to be degrees). For SOHO, this target phase 
had the additional attraction that it yielded a Lissajous 
that would have a halo-like trace in the jc-plane, 
avoiding the SEZ for at least a revolution or two. For 
Case AY3, the transfer involves a full revolution return 
to the Earth-side xr-crossing, so the phase angles needed 
were fay = 180 and y/ z = +90. Also different for Case 
AY3 is the introduction of a r-axis impulse located at 
the -^-extremum, which is mid-way between the 
departure and insertion bums. The fourth case, AY4, 
involves a !4 revolution transfer to the Sun-side to 
achieve an even smaller Lissajous, i.e.. Ay = 300,000 
km, with insertion state specified by A z = 120,000 km, 
fay = 0, and y/ z = -90. 

The four targeting schemes — a unique scheme for each 
of the four Ay cases — are given in Table 1. In the 
second column of Table 1, D = Departure, I = Insertion, 
and Z = r-axis bum. The Cartesian velocity components 
are denoted Vx, Vy, and Vz. To interpret the table, we 


Table 1. Ay Reduction Targeting Schemes 


Case 

AV 

Impulse (AV) 
Variables 
(components) 

Goal 

Variables 
at xz-plane 

Goal 

xz-plane at 
Impulse plus 

AY1 

D 

x,y,z 

Vx, Vy, Vz 

Zi rev 

I 

y 

o 

II 

£ 

1 rev 

AY2 

D 

x, y 

x,y 

V 2 rev 

I 

x, y, z 

Vx, Vy, Vz 

n/a 

AY3 

D 

*>y 

X,y 

V 2 rev 

Z 

z 

z 

Va rev 

I 

x,y,z 

Vx, Vy, Vz 

n/a 

AY4 

D 

x,y 

X,y 

V 2 rev 

I 

x,y 

Vx, Vy, 

n/a 


n/a = not applicable 


take Case AY3 as an example. Recall that the departure 
burn is located at the Earth-side xr-plane and is targeted 





Fig. 7. Case AZ4 RLP xr-projection 



urgent operations problem. The goal was to develop 
practical transfer techniques using existing software and 
methodologies. Since the transfer arcs developed have 
flight times between xr-plane crossings comparable to 
that of the original halo (within 3.3 percent) and achieve 
virtually tangential approaches to the target LPOs, they 
at least must be considered reasonably efficient. 

An important point relevant to practical operations is 
that SOHO — as well as most spacecraft — would need to 
perform these impulses in components. For example, 
for SOHO the Case AY2 cost rises from 55.9 to 71.2 
m/sec when performed in components. Likewise for 
Case AZ4, the cost increases from 106.2 to 1 15.1 m/sec. 

In the end, SOHO’s orbit was not changed thanks to a 
less risky alternative solution to the contingency having 
been found. Yet the LPO transfer techniques described 
here have potential applicability beyond just SOHO to 
virtually any LPO mission. Though the cases examined 
were limited in number and scope, it seems likely that 
these transfer techniques could be just as successful for 
modifying the y- and r-amplitudes in either direction 
and achieving any set of desired amplitudes. 
Applicability to L2 and lunar LPOs is also expected. 
However, these suppositions beg further research. 
Finally, in cases where A z amplification is desired, the 
suitability of a Lissajous versus a periodic halo should 
be investigated, as the delta-V costs could be less. 


Fig. 8. Case AZ4 RLP xy-projection 
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a periodic halo) are likely best for general applicability. 
For A y modification, however, if somewhat looser r-axis 
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that larger issue was not explored. The techniques were 
developed during a short-window effort to address an 
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